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ABSTRACT 



In 2007, R Coronae Borealis (R CrB) went into an historically deep and long 
decline. In this state, the dust acts like a natural coronagraph at visible wave- 
lengths, allowing faint nebulosity around the star to be seen. Imaging has been 
obtained from 0.5 to 500 /im with Gemini/GMOS, HST/WFPG2, Spitzer/MIPS, 
and Herschel/SPIRE. Several of the structures around R CrB are cometary glob- 
ules caused by wind from the star streaming past dense blobs. The estimated 
dust mass of the knots is consistent with their being responsible for the R CrB 
declines if they form along the line of sight to the star. In addition, there is a 
large diffuse shell extending up to 4 pc away from the star containing cool 25 
K dust that is detected all the way out to 500 /im. The SED of R CrB can 
be well fit by a 150 AU disk surrounded by a very large diffuse envelope which 
corresponds to the size of the observed nebulosity. The total masses of the disk 
and envelope are 10~^ and 2 M0, respectively, assuming a gas-to-dust ratio of 
100. The evidence pointing toward a white-dwarf merger or a final-helium-shell 
flash origin for R CrB is contradictory. The shell and the cometary knots are 
consistent with a fossil planetary nebula. Along with the fact that R CrB shows 
significant Lithium in its atmosphere, this supports the final-helium-shell flash. 
However, the relatively high inferred mass of R CrB and its high fluorine abun- 
dance support a white-dwarf merger. 



Introduction 



R Coronae Borealis (R CrB) was one of the first variable stars identified and its bright- 



ness variations have been monitored since its discovery over 200 years ago (Pigott & Engle- 



field 1797). The RCB stars form a small group of carbon-rich supergiants which are defined 



by extreme hydrogen deficiency, and unusual variability, characterized by large declines of 8 
mag or more due to the formation of carbon dust at irregular intervals ( Clayton]|1996 ). Two 
evolutionary scenarios have been suggested, a merger of a double-degenerate white dwarf 



merger, or a final helium shell flash (e.g., Webbink 1984 Clayton et al. 2007). The unex- 
pected discovery that some of the RCB stars have isotopic abundances of ^^O which are up 



^Based in part on observations with the NASA/ESA Hubble Space Telescope obtained at the Space 
Telescope Science Institute, and from the data archive at STScI, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555 
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to 500 times those seen in other stars has increased the hkehhood that these rare stars are 
the resuh of a white-dwarf merger dClayton et aL][2005| [20071 |Garcia-Hernandez ^t aL][20lo| . 

RGB stars have long been known to show IR excesses due to warm circumstellar dust 
(Clayton 1996 Feast et al. 1997, and references therein). Spectropolarimetry of R CrB 
obtained in a deep decline showed optical depth effects with wavelength that are consistent 
with a bipolar geometry with a thick disk or torus obscuring the star and additional more 
diffuse dust around the poles ( Stanford et al.||1988 Clayton et al.||1997 ). Resolved nebulosity 
has been detected from the UV to the IR around several RCB stars, including R CrB itself, 
ranging in size from a few arcseconds (V854 Cen, UW Cen) to arcminutes (R CrB, RY Sgr, 
V CrA, UW Cen) ( |Walker| [l985| |Rao fc Nand^[l986| [Walker] [l986| [PoUacco et aL][l99l 



Clayton et al. [1999 Clayton & Ayres 2001), possibly linking these stars with final helium 



shell flash stars and the central stars of planetary nebulae (PNe). In particular, R CrB has a 
far-IR shell which is ~20' in diameter, corresponding to 8 pc at its assumed distance ( Gillett 



et al.|1986 ). More recent observations have used interferometry and adaptive optics to obtain 



information about the distribution of dust around RCB stars on small spatial scales of tens 
of AU dOhnaka et aL][2003[ [de Laverny fc Mekarnia|[2004l [Leao et al.|[2007[ [Bright et al.][20llj ) 

In 2007 July, R CrB went into a historically deep and long decline reaching V=15. In 
this state, R CrB acts as a natural coronagraph, as the thick dust cloud in front of the 
star eclipses its photosphere but not its circumstellar material. In this paper, we present 
new Germm/GMOS, HST/WFPG2, Spitzer/MIPS, and Herschel/ SPIRE images of R CrB 
obtained just before or during this decline, giving wavelength coverage from 0.5 to 500 /im. 
We will investigate the morphology and nature of the dust clouds seen near R CrB as well 
as its large far-IR shell. 



2. Observations and Data Reduction 

The new observations obtained for this study are summarized in Table 1. Figure [T] shows 
the recent AAVSO lightcurve of R CrB with the epochs of the various observations marked. 

R CrB was observed with Gemini-S/GMOS on 2009 March 23. Images were obtained 
in the g', r', i', z', and CaT filters. We obtained 8 integrations of 30.5 s in each of the g' and 
r' filters for a total of 244s, and 14 integrations of 30.5 s each in the i', z', and CaT filters for 
a total of 427s. We used the 2x2 binning and slow readout mode. The field is 6'x 6'. The 
central part of the field is shown in Figure [2j All images have been geometrically registered 
to within 0.1 pix rms. A point spread function (PSF) centered on R CrB has been subtracted 
where the sky brightness has been ramped down to zero as one goes out 27 pixels in radius. 
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The plate scale is QfHAG per pixel. The seeing was ^&!8. These PSF-subtracted images are 
shown in Figure |3| 

R CrB was observed with HST/WFPC2 on 2009 April 27. Images were obtained with 
the F555W and F814W filters with total integration times of 487s and 496s, respectively. The 
individual images were stacked and drizzled. The images were then registered and a PSF was 
subtracted. These images are shown in Figure |4j WFPC2 images of R CrB taken on 1996 



July 1 (JD 2450265) with the F469N filter were downloaded from the MAST archive (Ueta 



eF"aL][2000l ). The exposures were 4 x 20s and 2 x 400s. When these images were obtained 
R CrB was at a visual magnitude of ~8.1 (AAVSO) while recovering toward maximum light 
after a deep decline. The individual images were stacked and drizzled. 



A UBVRI photometric sequence has been established for the R CrB field (Landolt & 



ae5^||2010j ). A number of stars (Nos. 13, 14, 19, 20, 22, 23, 31, 34, 35, 37, 39) from that 
sequence are also in the Gemini/GMOS field. Two of the stars (Nos. 20, 22) lie in the 
HST/WFPG2 field but none are on the PC chip. In addition, a strip from the Sloan Digital 
Sky Survey (SDSS) telescope crosses R CrB. Several of the photometric sequence stars (Nos. 
13, 14, 19, 20, 23, 35, 37, 39) from [Landolt fc Clem| ( [2010l ) also have ugriz photometrjQ and 
lie in the CMOS field. 

Spitzer/MIPS scan maps were obtained of R CrB at 24, 70, and 160 /im on 2007 April 
14/15 as part of program 30029. We mapped a 20' wide area using 15 legs at the medium scan 
rate with 80" (1/4 array) cross-scan offsets. The total integration time was 15,536s. R CrB 
was observed twice to allow for good Germanium (70 and 160 /im) transient and asteroid 
removal. The second epoch maps allow for the clean removal of the Germanium responsivity 
drifts which are especially important for mapping extended emission around R CrB. The 



scan maps were reduced and mosaiced using the MIPS Data Analysis Tool (Gordon et al. 



2005) supplemented with custom reduction scripts written specifically to improve the MIPS 



reductions of extended sources. The custom reduction scripts include extra stpdf beyond 
that of the MIPS DAT. At 24 /im, the extra stpdf include readout offset correction, array 
averaged background subtraction (using a low order polynomial fit to each leg, with the 
region including the Galaxy excluded from this fit), and exclusion of the 1st five images in 
each scan leg due to boost frame transients. At 70 and 160 /im, the extra processing step 
is a pixel dependent background subtraction for each map to remove residual detector drifts 
and background cirrus and zodiacal emissions. The 24 /im image was saturated in the core 
and was therefore excluded from the photometric analysis. 



Herschel/SPIRE (Pilbratt et al. 2010 Griffin et al. 2010) images at 250, 350 and 500 



^http; / / www.sdss.org/ 
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/iin were obtained on 2009 December 27 with a total integration time 4408s as part of the 
Herschel key programme, "MESS - Mass-loss of Evolved StarS." The reduction of the SPIRE 



data is described in Groenewegen et al. (2011). The data were processed using the standard 



Griffin et al. 


2008, 


2010 


Dowell et al. 


2010 



The procedure for calibrating the data is described in Swinyard et al. (2010). 



The IRAS imaging data have been reprocessed using a better zodiacal light subtraction 
and better destriping (Miville-Deschenes & Lagache 2005). The co-added 100 /xm IRIS 



images have been downloaded for use in this study. The IRAS, Spitzer and Herschel images 
are shown in Figure |5j The new images have significantly better resolution than the IRAS 100 
fim images. The pixel sizes are IRAS/100 (1^5), MIPS/70 (5"), MIPS/160 (16"), SPIRE/250 
(6"), SPIRE/350 (10"). SPIRE/500 (14"). 



Discussion 



R CrB lies at high Galactic latitude (b =-|-51°) and so the foreground extinction is 



quite small, E(B-V)~0.035 mag (Schlegel et al. 1998). At maximum light, R CrB is V=5.8 



mag and B-V =0.6 mag ( Lawson et al.||1990[ ). The absolute magnitude of R CrB is estimated 
to be My =-5 mag based on the absolute magnitude/effective temperature relationship found 
for the Large Magellanic Cloud RGB stars ( [Alcock et aL]|200H|Tisserand et al.||2009| ). So, for 
the analysis in this paper, we adopt a distance of 1.4 kpc. The most recent decline of R CrB 
began in 2007 July (~JD 2454290) and the star began returning toward maximum light in 



the spring of 2011. Landolt & Clem (2010) point out that the recent decline of R CrB is 



among the longest ever recorded. This decline was deeper than most, bottoming out at V 
~15.2 mag. 



3.1. Background Objects 

As shown in Figure [l| The Gemini/GMOS and HST/WFPG2 imaging was obtained 
when the star was at V~15 mag or about 9 magnitudes below maximum light. The extreme 
faintness of the star revealed for the first time that R CrB lies in front of a cluster of galaxies, 
shown in Figure |2} 

A color magnitude diagram (CMD) was constructed using the Gemini/GMOS images. 
The CMOS g-band image was registered to the r-band image. Then SExtractor was used 
to obtain photometry in dual image mode, with the r-band image used for selecting objects. 
For the g-r colors in the CMD, photometry was obtained through 2" diameter apertures. 
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The GMOS g — r vs r diagram, shown in Figure [6| displays a fairly clear red sequence at 
g — r ~ 1.6 and r ~ 20 which corresponds to a cluster of galaxies at z ~ 0.3. The sizes of 
the galaxies in the GMOS imaging are appropriate for this redshift. 

R CrB also lies in a SDSS strip obtained in 2005 when the star was bright. SDSS did not 
obtain spectra of any of the galaxies in Figure [2| but it did get spectra of about 25 galaxies 
within ~10' of R CrB. Of these 14 had redshifts in the range 0.07-0.09, and 6 had redshifts 
in the range 0.16-0.18. The SDSS spectroscopy is unable to reach normal cluster galaxies at 
z ~ 0.3 and the SDSS redshifts are likely to be those of field galaxies in the foreground of 
the cluster. 



Landolt & Clem (2010) were the first to point out that there is a faint red star ~3" west 
of R CrB. They find that the star has a magnitude of r=20.78±0.06 and 5-1/ = 0.99±0.08. 
It is also clearly visible in the HST, and Gemini images taken during the R CrB minimum in 
Figures |3] and |4j If the star is a physical companion, it would be at a projected distance of 
about 4000 AU from R CrB, and its absolute magnitude would be ~+10 mag. This would 
imply that the star is a main sequence M star, however the observed B — V colors are not 
consistent with such a star. On the other hand, if the faint companion to R CrB is a K2 V 



star, which is more consistent with its colors {B — V = 0.91 and My=+6.4 mag ( Cox||2000 )) 



then it is at a distance of ~7 kpc. In short, the companion star is likely to be a background 
K dwarf. The companion star to R CrB is not visible in the HST/WFPC2 F469N images 
from 1996. The WFPC2 ETC, assuming a 20th magnitude KO star, estimates that the ratio 
of the count rates in the three filters will be F469N/F555W/F814W = 0.04/5.7/6.7 electrons 
pixel"^ s^^. So it is likely that the star was too faint to be seen in the F469N images. 



3.2. Cometary Knots Close to R CrB 

The Gemini and HST fields centered on R CrB are shown in Figures [3] and |4j A stellar 
PSF has been subtracted from the position of R CrB for each image. The Gemini and HST 
images were obtained only a month apart and the same structures can be seen in both. The 
structures are labeled in Figure |4| showing real nebulosity as well as instrumental artifacts 
such as diffraction spikes and charge transfer effects. The real structures appear in multiple 
filters in both GMOS and WFPC2 images so are unlikely to be artifacts. 

R CrB was 9 magnitudes fainter than at maximum light during the decline. The dust 
lying in front of the star is extremely optically thick so at this time we are seeing essentially 
zero directly transmitted light. When the star itself is completely obscured, all that remains 
is a fairly constant scattered-light residual that is likely to be from photons that have escaped 



around the edge of the obscuring cloud by scattering off dust lying beside or behind the star. 
This scattered light may be responsible for the typical flat-bottomed declines often seen in 
RGB star lightcurves. When the stellar photosphere is obscured in deep declines, RGB stars 
show emission lines in their spectra. Spectra were obtained of R GrB in 2008 when it was 
deep in the present decline showing emission lines including the Na I D lines and the IR Ga 



II triplet (Rao & Lambert 2010). These spectra show that the emission lines account for no 



more than 5-10% of the flux in the V band. R GrB, itself, appears to be marginally resolved 
in the WFPG2 images. There is a diffuse, nearly circular "halo" around the star. This may 



be the dust near the star that is scattering the light around the obscuring cloud. Jeffers et 



al. (2011) also suggest the presence of this halo based on the brightness of R GrB at the 



bottom of its decline. The blobs visible in Figures [3] and |4] are also seen in scattered light 
from R GrB. 

Several of the structures around R GrB appear as jet-like blobs pointing radially away 
from the star. This could either be a lighting effect created by the central star shining 
through holes in the dust clouds and illuminating different portions of the circumstellar 
dust shell, or the blobs could be "cometary" knots. The RGB star, UW Gen, shows both 



illumination structures and cometary structures (jGlayton et al. 1999| Glayton et al. 2011, 
in preparation). 



Figure [7| shows a comparison of the blobs seen in R GrB and cometary 
globules seen in UW Gen. In particular, the globule to the SE of R GrB shows the typical 



"tadpole" shape similar to those seen in the Helix nebula (Matsuura et al. 2009). It also 



shows a sharp edge on the side toward the star with the brightness fading away from the 
star, as shown in Figure |8} The largest knot (#1 in Figure [T]) is ~1000 AU across. This 
is larger than the knots in the Helix which are 100-300 AU in size with masses of 1O~^M0 



(O'Dell & Handron 1996). The simplest explanation for these cometary knots is that they 



are caused by a stellar wind streaming past a density enhancement moving slowly away from 



the star ( Garcia-Segura et al. 2006). There is evidence for both low (10-20 km s ^) and 



high (200-400 km s ^) velocity gas streaming away from RGB stars (Glayton et al. 1992 



2003 Garcia- Hernandez et al. 2011). Simulations show that when a slow wind is overtaken 



by a fast wind, the gas-gas interaction could create knot- like shapes (Pittard et al. 2005). 



Recent investigations using a hydrodynamics code including dust shows that small dust 



grains (<0.045 jdm) will follow the motion of the gas closely (van Marie et al. 2011). The 
UW Gen nebula, about 15" across, resembles a fossil PN shell where the gas is neutral and 



the shell is seen in light scattered from dust (Glayton et al. 1999). UW Gen is at a distance 
of ~5.5 kpc so if it were at the distance of R GrB, its nebula would be about 1' in size. As 
described below, the diffuse nebulosity around R GrB extends at least 10' from the star but 
no blobs or density concentrations are seen more than a few arcseconds from the star. 



Jeffers et al. 



(2011) report imaging polarimetry of the cometary knot SE of R GrB (^^1 



- 8 - 



in Figure [7]). Combining those data with the WFPC2 images presented here, they infer that 
the knot is at a scattering angle of 105° and has a dust mass of ~10-6 M0. It hes about 
1" from R CrB which is about 1400 AU from the star. At velocities of 20 and 200 km s~^, 
the age of the blob would be 300 or 30 yr, respectively. 

To estimate the mass of dust present in the knots, we make the following simplifying 
assumptions: (1) the illuminating flux from R CrB is unobscured (i.e. equal to the maximum 
values given in Table 2), (2) the dust in the knot is optically thin, such that all grains are 
illuminated and no grains are shadowed, and (3) that the grains can be described by the 



neutral PAH and graphitic solid properties from iLi & Draine (2001). Following Chevalier 



(1986), the flux scattered off one dust grain of radius a at a distance r from a source with 



flux F(A) is 

dFscaiX, r, a) = — (1, 



where Qsca is the scattering efficiency, $ is the Henyey & Greenstein (1941) phase function, 
and 6 is the scattering angle. If one assumes the dust grains have uniform size, the mass of 
grains {Mg^) in a knot that produces a total scattered flux (Fgca) is 

3FiX)QUKaMX,a,e) ^ ^ 

where p = 2.24 g cm~^ is the density of carbonaceous dust. If one instead adopts an MRN 
size distribution ( Mathis et al.||l977 ) such that the number density of grains Ugr oc a~^^^ then 



F(A)/0.2247ra-i-5Q,,„(A,a)<l>(A,a,e)da ^ ' 



where the numeric constant is the normalization factor from Weingartner & Draine (2001), 
and assuming a gas-to-dust ratio of 100. Finally, if the dust is assumed to lie in a spherical 
shell of inner and outer radii Rm and Rout then Mgj. must be integrated over the volume of 
the shell. Since 



$(A,a,^)sin^rf0 = 2 (4) 



the total grain mass for a single-sized model equals 

9F{X)Qsca{X, a) Rout — Rin 



Of the knots identified in Fig. 4, numbers 4 and 6 are coincident with diffraction spikes, 
and we refrain from estimating their dust masses. For the others, the total fiux from each 
knot was measured through manually-sized elliptical apertures. These values are reported in 
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Table 3. For each knot, models with single grain sizes of 0.001, 0.01, 0.1 and 1.0 /im, as well 
as full MRN distributions, were tested for scattering angles from 10°to 170° (in increments of 
10°) to find which combinations of grain sizes and scattering angles yielded equivalent dust 
masses for both the F555W and F814W photometry. The best-fitting models, including the 
grain sizes, scattering angle, resulting clump distance from R CrB and estimated dust mass, 
are also listed in Table 3. The model dust masses are 10~^° to 10~^ Mq. Note that if a 
model is not listed, it failed to reproduce the photometry in both filters. 

Knot 3 has very low flux in the F814W filter, resulting in a large uncertainty with the 
flux between 4 and 6 xlO~^® erg cm~^ s~^ A~^. For this reason, results are shown for the 
upper and lower flux range. Since small grains scatter isotropically, there were no preferrable 
scattering angles for this knot using the smaller single grain-size models, making a distance 
estimate from R CrB impossible. Owing to the very bright F814W flux from knot 5, we find 
there is no dust-scattering model that adequately reproduces the fluxes in both filters. 

To estimate the dust mass producing the halo (i.e., the slightly- resolved profile around 
R CrB), we assume that all of its flux is light scattered by a thin circumstellar shell. We 
find that a limb-brightened shell of optically-thin dust can be up to 0'.'3 in radius and still 
produce a profile consistent with that of R CrB in the 2009 HST images. The flux of the 
halo, as well as the best-fitting dust model that produces this scattered light, is listed at the 
bottom of Table 3. 

Note that our first modeling assumption implies that the masses derived above are lower 
limits, since more dust is required if the incident flux is lower than that of R CrB at maximum 
light. As a quick check of the second assumption, we assume the knots are ellipsoids with 
the same axes as those used to measure the fluxes. In knot 1, there are 8 x 10^^ grains of 
radius 0.1 fim, corresponding to a number density of n^^ = 8 x 10~^° cm~^. The optical 
depth T = ngrL{Qsca + Qabs)T^o? where L is the line-of-sight distance through the dust and 
Qahs is the absorption efficiency. The proposed ellipsoid has a depth around 10^^ cm, yielding 
an approximate l^-band optical depth of 8 x 10~^, which is sufficiently thin to warrant the 
assumption that all grains are illuminated by the star. See §3.3 for a justification of the last 
assumption, namely, that the dust is purely carbonaceous. 



3.3. The Large Far-IR Dust Shell 



R CrB has an IRAS sh.e\\ with a radius of 10', 4 pc at a distance of 1.4 kpc (Gillett et al. 
1986). The new MIPS images at 70 and 160 yum, shown in Figure [s] along with the IRAS 100 
/xm image, confirm the presence of a large dust shell. It is also visible in the SPIRE images 
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at 250 and 350 fim. It is not clearly visible at 500 fim. The new MIPS and SPIRE images 
have significantly better resolution than IRAS. These new images imply that the shell is 
spherical in shape and that the component to the northeast of the nebula in the IRAS 100 
/xm image is detached from the R CrB shell in the MIPS and SPIRE images. The SPIRE 
images, in particular, are full of point sources, which are background galaxies. To determine 
how much of the nebulosity seen in the IRAS, MIPS and SPIRE images might be due to 
blended emission from background galaxies, we used Sextractor to identify and remove point 
sources from the SPIRE 250 fim image. The result is shown in Figure [9] This clearly shows 
that the dust nebulosity around R CrB is not due to a concentration of galaxies in a cluster 
behind the star. There is significant diffuse far-IR emission, centered on R CrB. To get a 
clearer idea of the extent of the nebulosity at different IR wavelengths, the average brightness 



in 1' concentric rings from R CrB was calculated and plotted in Figure 10 These new data 



show that the shell is possibly even larger than that seen by IRAS. The MIPS 70 /im image 
also has a central elongated structure running almost east-west that has a diameter of ~3'. 

At the assumed distance of R CrB, the very extended IR shell extends about 8 x 10^ 
AU from the star. The dust mass and total mass of this shell have previously been estimated 
to be ~10-3-10-2 M0 and ~0.3 - 3 Mq, respectively (|Gillett et aLl|1986[). They estimate 



that the dust temperature in the shell is 25-30 K. Rao & Nandy (1986) find a hot component 



close to the star with a dust mass of 5 x 10 ^ Mq and a cool component with a dust mass 
of 10-2 - 10-3 M. 



10. 



Photometry available from the literature for R CrB has been combined with new pho- 
tometry using the MIPS and SPIRE images. The data, summarized in Table 2 and plotted 



in Figure 11 consist of visible photometry at maximum light, JHKL, as well as IRAS pho- 



tometry (Walker 1985; Rao fc Nandy 1986). AKARI/IRC and FIS photometry are also 



available (Murakami et al. 2007; Ishihara et al. 2010). Fluxes were extracted from SPIRE 
(250, 350, 500 fim) using Starfinder (Diolaiti et al. 2000). Starfinder produced PSF fit- 
ted photometry using PSFs which were downloaded from the Herschel Science Center. No 
aperture corrections were applied since the PSFs used were of sufficiently large diameter. 
Aperture photometry was also performed as a test of the Starfinder results. The photomet- 
ric uncertainties quoted in Table 2 for SPIRE are 15%, or the Starfinder uncertainties plus 



calibration errors of 5% whichever is higher (Swinyard et al. 2010). 



We have modeled the SED of R CrB using a Monte Carlo radiative transfer code which 
includes nonisotropic scattering, polarization, and thermal emission from dust in a spherical- 
polar grid dWhitney et"aL]|2003a||b| [Robitaille et~aL]|2006D . All of the models were done with 
amorphous carbon dust and an MRN size distribution ( Mathis et al.|1977 ). We have assumed 
a mass gas-to-dust ratio of 100. The best fit to the R CrB SED is shown in Figure [TTl A 
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good fit was obtained to both the hotter dust centered on the star, as well as the cooler, 
large diffuse nebula. Heating by the interstellar radiation field is included. The best fitting 
model is a small disk surrounded by a large diffuse envelope. The disk has inner and outer 
radii of 42 and 150 AU, and the envelope has inner and outer radii of 1.3 and 8.8 x 10^ AU. 
The envelope has a bipolar cavity carved out with an opening angle of 30°. The Ay through 
the envelope is quite low, less than 0.01 mag over a wide range of viewing angles (0°-75°). 
The disk dust mass is 3.5 x 10~^ Mq and a total mass of 3.5 x 10"*^ Mq. Assuming a 
density profile of r^^ which is appropriate for dust accelerated away from the star radially 
by radiation pressure, the dust mass of the outer envelope is 2 x 10~^ Mq and the total 
mass (dust + gas) is 2 Mq. The gas may be dragged along with the dust or it may become 
decoupled. 



3.4. The Evolution of R CrB and its Dust 



The generally accepted model for dust formation in the RGB stars, the condensation of 



carbon dust along the line of sight, was suggested over 70 years ago (Loreta 1935 O'Keefe 



1939 Clayton 1996). The evolution of RGB star spectra and hghtcurves during declines is 



consistent with dust that forms close to the stellar atmosphere and then is accelerated to 



hundreds of km s ^ by radiation pressure (Glayton et al. 1992; Whitney et al. 1993). If the 
dust forms near the star, it will accelerate very quickly, possibly reaching 200 km s~^ 



m 



30 days. Other evidence supports a much slower expansion (e.g., Garcia-Hernandez et al. 



2011, and references therein). The observed timescales for RGB dust formation fit in well 



with those calculated by carbon chemistry models ( 


Feast 


1986; 


Woitke et al. 


1996 


Rao & 



Lambert 2010). 



The dust forming around R GrB is not in a complete shell but rather in small "puffs" 
(Glayton 1996). When a puff forms along the line of sight to the star we see a decline. 



Estimates of the covering factor of the clouds around RGB stars during declines, from both 



extinction studies and IR re-emission of stellar radiation, indicate that f < 0.5 (Feast et al. 



1997 Glayton et al. 1999 Hecht et al. 1984 1998). The typical puff causing a decline is 



thought to have a dust mass of ~10-s Mq ( |Feast||1986| [Glayton et aLl|1992[ ). There is strong 
evidence that the puffs form in or near the surface of the RGB star due to density and 



temperature perturbations caused by stellar pulsations (Woitke et al. 1996 Grause et al. 



2007). The RGB stars show regular or semi-regular pulsation periods in the 40-100 d range 



( Lawson et al.|[l990 ). R GrB, itself, does not have one regular period but has shown periods 



of 40 and 51 d (Fernie & Lawson 1993). 



For this recent deep decline, with Av=9 mag, assuming the dust forms at 2 Rj, (Rj,=85 
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Rq), and a puff subtends a fractional solid angle of 0.05, then the dust mass will be ~10~^ 
Mq. The puff would be accelerated quickly by radiation pressure and would dissipate rapidly 
so dust must be formed continually by R CrB to maintain itself in a deep decline for 4 years. 
If a puff is formed during each pulsation period of ~40 d, R CrB would be forming about 



10"^ Mq of dust per year. Assuming a gas-to-dust ratio of 100 rtWhittet 2003L the total 
mass loss per year is 10~^ Mq. The estimated masses of the puffs causing the declines and 
the cometary knots is similar. The puffs cause declines when they form directly in our line 
of sight and may be seen as cometary knots when they form to the side of or behind R CrB. 

Little is known about the lifetime of the RCB phase. We have a lower limit from the fact 



that R CrB itself was discovered to have large brightness variations 200 years ago (Pigott & 



Englefield 1797). Assuming a velocity of 20 km s~^ then the large diffuse dust shell around 
R CrB would take 10^ yr to form. If the mass-loss was more like the high-velocity winds 
seen in R CrB today (~200 km s~^) then the shell would be about an order of magnitude 



younger ( Clayton et al.|[2"003 ). If R CrB is the result of a final- helium-shell flash rather than 
a white-dwarf merger, then the size and timescales would be consistent with the nebulosity, 
now seen in far-IR emission, being a fossil PN shell. The nebulosity including cometary 
knots, seen around R CrB and UW Cen, shown in Figure [7} looks very much like a PN shell 



(Clayton et al. 1999, Clayton et al. 2011, in preparation) 



Planetary Nebulae are a very late phase in the evolution of low and intermediate mass 
[1-8 Mq) stars. The PN, NGC 6302 is estimated to have an envelope with a total mass of 



4.7 Mq and a central star with a progenitor mass of 5.5 Mq (Wright et al. 2011). So the 



estimate made here for the total envelope mass of R CrB is consistent with it being a fossil 
PN shell. Any gas lost during the white-dwarf merger would have far less mass. If the shell 
is an old PN shell then this would suggest that R CrB is the product of a final flash event 
rather than a white-dwarf merger. The estimated mass of the envelope is consistent with a 
mass-loss rate of 10^^ Mq for about 10^ yr. 

Five Galactic RCB stars, including R CrB, itself, show significant Lithium in their at- 



mospheres (Rao & Lambert 1996 Asplund et al. 2000 Kipper & Klochkova 2006). Renzini 



(1990) suggested that in a final flash the ingestion the burning of the H-rich envelope leads 



to Li-production through the Cameron- Fowler mechanism (Cameron & Fowler 1971). The 



abundance of Li in the atmosphere of the final-flash star, Sakurai's object, was actually 
observed to increase with time ( Asplund et al.|[T999 ). In a white-dwarf merger, the tempera- 
tures required for the production of ^^0 should destroy any Li present. So, the simultaneous 
enrichment of Li and ^^0 is not expected in the white-dwarf merger scenario. The abundance 
of ^^O cannot be directly measured in R CrB, but it is overabundant in fluorine, which does 
imply a high ^'^O abundance (Pandey et al. 2008). 
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About 10% of single stars will undergo a final-fiash event ( Iben et aL]|1996 ). About this 
percentage of RGB stars (R CrB, RY Sgr, V CrA, and UW Cen) show evidence of resolved 
fossil dust shells ( Walker|1994 ). R CrB is thought to be ~0.8-0.9 M0 from pulsation modeling 
(Saio 2008), and this mass agrees well with the predicted mass of the merger products of a 
CO- and a He- WD (Han 1998). On the other hand, final-flash stars, since they are single 
white dwarfs, should typically have masses of 0.55-0.6 Mq (Bergeron et al. 2007). 



4. Summary 

New images were obtained of R CrB with Germm/GMOS, HST/WFPG2, Spttzer/MIPS, 
and Herschel/SPIRE just before or during its recent very long and very deep decline. At 
visible wavelengths with Gemini and HST, the dust obscuring R CrB acts as a natural 
coronagraph, allowing material close to the star to be imaged for the first time. In particular, 
several cometary knots have been detected within a few thousand AU of the star. The puffs 
of dust that cause declines in R CrB and the cometary knots have similar inferred masses 
so they may be the same dust clouds seen in and out of the direct line of sight to the star. 
Also, the existence of a cluster of galaxies directly behind R CrB was revealed for the first 
time. The IR images confirm the existence of a very large diffuse shell (radius ~4 pc) around 
R CrB that contains cool dust (~25 K) which is detected out to 500 fim and has a total 
mass of ~2 M0. This places a possible upper limit on the lifetime of R CrB as an RCB of 
~10^ yr. In addition, there is no significant mass of cold dust around R CrB that was missed 
by IRAS. 

RCB stars are rare so they don't contribute significantly to the interstellar medium, 
but their continual dust formation makes them ideal laboratories for understanding dust 
formation in evolved stars. In particular, both the RCB stars and Type la supernovae may 



Wheeler 


1978 


Webbink 


1984; 



Kilic et al. 2010). So studies of the circumstellar environments of RCB stars may provide 
information on the material around Type la supernovae. It is possible that there are two 
avenues for producing RCB stars and that R CrB is a product of the final flash rather than 
the white-dwarf merger. But the evidence is mixed. The IR shell and the presence of lithium 
point toward a possible final flash, while the relatively high mass of R CrB and the presence 
of fluorine point to a WD merger. 
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Fig. 1. — AAVSO data for R CrB since 2004. Visual magnitudes are plotted as black dots. 
Johnson V data are plotted as blue open circles. The epochs at which the Spitzer, Gemini, 
HST and Herschel data were obtained are marked. 



- 20 - 




Fig. 2. — Gemini/ GMOS color composite image with g' (blue), r' (green), and z' (red), 
showing the large number of galaxies near R CrB which is the bright object in the upper 
left. The field is 2'x 2'. North is up and East to the left. 
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Fig. 3. — Gemini/GMOS g', r', i', z', and CaT images. A stellar PSF has been subtracted 
from the position of R CrB for each image. There is a faint star 3" west of R CrB and a 
galaxy about 15" northeast of R CrB. North is up and East to the left. 



- 22 - 



F555W I 

I F814W 



Ring 5 

i \ 



2 

6 



Diff. Spike 



Ciiarge Transfer^ 

\ ^ 
Diff. Spiite -; 



Fig. 4.— i/5r/WFPC2 F555W (left) and F814W (right) images with a PSF removed from 
the position of R CrB. White arrows indicate instrumental features and the companion star. 
Blue arrows indicate real nebular structures around R CrB. The angular distance between 
R CrB, which lies in the center of the images, and the faint star to the west is ~3". The 
cometary knots are numbered. 
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Fig. 5. — Far-IR images centered on R CrB showing the large shell with a radius of 10'. 
Upper: IRAS/IRIS 100 /im, Spitzer/MIPS 70, and 160 /im, Lower: Herschel/ SPIRE 250, 
350, and 500 /xm. The field is 35' x 35'. 
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Fig. 6. — The CMD based on the Gemini/GMOS images. There is a clear red sequence at 
g-r = 1.6 which suggests a cluster at z~0.3. The clustering at g-r~l is probably M stars 
and foreground galaxies. 




Fig. 7.— HST/WFPC2 F555W image (17" x 17") of R CrB (left) and WFPC2 606W image 
(12" X 12") of UW Cen (right). The inset shows the cometary globule (#1) to the lower left 
of R CrB. 
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Fig. 8. — A cut from the center of the PSF-subtracted through the large cometary knot 
(#1) in the HST/WFPC2 F555W image. The cometary globule lies at 1000-2000 AU from 
the star. 




Fig. 9. — Left: Herschel/SPIKE 250 fim image centered on R CrB. Right: Same image with 
points sources subtracted to show diffuse IR emission is present. The field is 34' x 34'. 
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Fig. 10. — Surface brightness of diffuse emission surrounding R CrB in one arcminute rings. 
The curves have been shifted vertically for clarity. 
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Fig. 11.— Maximum light SED for R CrB. Filled red squares: UBVRIJHKLM, open red 
circles: IRAS (12, 25, 60, 100 /im), open green squares: MIPS (70, 160 yum), filled green 
squares: SPIRE (250, 350, 500 /im), open red triangles: AKARI 9, 18, 90 yum), blue line: 
ISO spectrum. The dust associated with R CrB was clearly detected in all the SPIRE bands. 
Also plotted is the best fitting radiative transfer model to the SED (purple dashed line). 
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Table 1. New Observations 



Date 


Telescope 


Instr. 


Filters 


2007 April 14/15 


Spitzer 


MIPS 


24, 70, 160 fim 


2009 March 23 


Gemini-S 


CMOS 


g', r', i', z', CaT 


2009 April 29 


HST 


WFPC2 


F555W, F814W 


2009 December 27 


Herschel 


SPIRE 


250, 350, 500 fim 



Table 2. R CrB Photometry 



Band" 


Flux (Jy) 


^ (Jy) 


U 


4.53eH 


-00 


3.90C-02 


B 


1.16e+01 


1.07C-01 


V 


1.79eH 


-01 


1.66e-01 


Rc 


2.05eH 


-01 


1.90C-01 


Ic 


2.07cH 


-01 


1.91C-01 


J 


1.77e+01 


1.65e-01 


H 


1.43e+01 


l.OOe-01 


K 


1.43eH 


-01 


1.33e-01 


L 


2.56eH 


-01 


2.36e-01 


AKARI/9 


5.30eH 


-01 


2.44e+00 


AKARI/18 


2.15eH 


-01 


2.90e-02 


IRAS/ 12 


3.89eH 


-01 


1.55e+00 


IRAS/25 


1.71eH 


-01 


6.84e-01 


IRAS/60 


3.94eH 


-00 


3.15e-01 


MIPS/70 


2.03eH 


-00 


3.4E-02 


AKARI/90 


1.49eH 


-00 


1.14e-01 


IRAS/100 


2.00eH 


-00 


1.60e-01 


MIPS160 


2.97E- 


■01 


9.36E-03 


SPIRE/250 


7.81E- 


■02 


1.17E-02 


SPIRE/350 


3.40E- 


■02 


5.10E-03 


SPIRE/500 


1.25E- 


■02 


4.34E-03 



''Maximum light UBVRI photometry from |Cottrell et al.| 
l |1990^ . JHKL photometry from [Feast et al.| i |1997f 
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Table 3. Cometary Knot Models 



Knot 


-^F555W 


-Pf814W 


Grain 


e 


r 


Mgr 




10~^^ erg cm 




Model 





AU 


Mo 


1 


723 


261 


0.1 yum 


no 


1500 


3.7 X 10-9 








MRN 


35 


2500 


1.6 X 10-^ 


2 


45 


12 


0.1 /im 


70 


4700 


1.2 X 10-9 


3 


32 


4 


0.01 urn 


varies 


N/A 


(1-3) X 10-^ 








0.001 /xm 


varies 


N/A 


(1-5) X 10-^ 






6 


0.1 fim 


50 


3000 


2.5 X 10-1° 


Halo 


15600 


10600 


1.0 iim 


N/A 


<400 


2 - 4 X 10-s 



